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Abstract 

The evaluation of Yukawa-enhanced two-loop contributions to the 
MSSM Higgs-boson mass is considered. We prove the common assumption 
that regularization by dimensional reduction preserves supersymmetry at 
the required level. Thus generating counterterms by multiplicative renor- 
malization is correct. Technically, we identify a suitable Slavnov- Taylor 
identity, use a recently developed method to evaluate it at the two-loop 
level, and show that it is valid in dimensional reduction. 



1 Introduction 

The prediction of the mass of the lightest Higgs boson is one of the most striking 
features of the Minimal Supersymmetric Standard Model (MSSM). At tree level, 
the mass of the lightest MSSM Higgs boson is bound to be smaller than the 
mass of the Z-boson. At higher orders, Mh is a function of all MSSM param- 
eters, but there still is a very restrictive upper bound as long as the masses of 
supersymmetric particles are not much higher than 1 TeV. 
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The LEP exclusion bound of 114.4 GeV [1] for the mass of a standard model 
(SM)-like Higgs boson allows to derive stringent constraints on the MSSM pa- 
rameter space. In the future, the direct measurement of the Higgs-boson mass 
will lead to very small experimental uncertainties of only 200 MeV at the LHC [2] 
and 50 MeV at the ILC [3]. Comparing the measurement of M h to its predicted 
value amounts to sensitive tests of the MSSM and allows the indirect determi- 
nation of further MSSM parameters. It is therefore important to know as a 
function of the MSSM parameters as precisely as possible. 

A lot of effort has been put into the evaluation of (see Ref. [4] for a 
review). So far, the one-loop contributions are completely known. The most 
important two-loop corrections are the ones enhanced by the Yukawa couplings 
of the top and bottom quark a tt b = Utb/i^ 71 )- The Yukawa-enhanced contribu- 
tions of 0{a t a s ) [5-13], 0{a 2 t ) [5,14,15], 0{a b a s ) [16,17], 0(a t a b ), 0{a 2 b ) [18] 
are known in the on-shell renormalization scheme. All these Yukawa-enhanced 
contributions can be evaluated in the gauge- less limit, where the electroweak 
gauge couplings and M w>z go to zero with fixed ratio M w /M z and fixed Higgs 
vacuum expectation values. Further electroweak two-loop contributions that go 
beyond the gauge-less limit are known in the DR scheme [19,20]. The remain- 
ing intrinsic error of the Yukawa-enhanced corrections has been estimated to 
5M™ tT = 3 GeV [21]. 

There might be, however, an additional Yukawa-enhanced contribution from 
non-supersymmetric counterterms at the two-loop level. So far, such terms have 
been assumed to vanish, but without an explicit investigation. In this paper we 
close this gap by a detailed evaluation of the counterterm structure at two-loop 
order. 

In all the references listed above, regularization by dimensional reduction 
(DRED) [22] has been employed and it has been assumed that all symmetries 
of the MSSM are preserved by this method. According to this assumption, the 
structure of the counterterms is symmetric and generated by multiplicative renor- 
malization. However, although dimensional reduction has been checked to pre- 
serve supersymmetry in many individual cases, there is no general proof. 

The traditional checks apply to one- loop self energies [23], one-loop on-shell 
[24] and off-shell [25, 26] three-point functions and further one-loop relations in- 
volving soft and spontaneous symmetry breaking [27]. Further checks concern 
the renormalization group (3 functions at the two- loop level [28], corresponding 
to divergences of two-loop diagrams. Most recently, supersymmetry identities for 
one-loop four-point functions and two-loop two-point functions (without spon- 
taneous symmetry breaking) have been checked [29]. In all cases, dimensional 
reduction was in agreement with the required supersymmetry relations. These 
checks are indeed sufficient to prove that multiplicative renormalization is correct 
for the one-loop counterterms of the Higgs, quark and squark sectors entering the 
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two-loop calculation of Mh at the level of one-loop subrenormalization. 

These results, however, do not constitute a proof that dimensional reduction 
preserves supersymmetry at the two-loop level. In particular, they cannot exclude 
that supersymmetry is broken by a finite amount, which could be relevant for 
precision calculations of M h . 

In the case of supersymmetry-breaking by the method of regularization, 
extra symmetry-restoring counterterms have to be added as discussed e.g. in 
Refs. [25-27]. These extra counterterms, which are by themselves not supersym- 
metric, are required to restore the original symmetry at the quantum level at the 
considered order. They are finite since all divergences are cancelled by symmetric 
counterterms derived from multiplicative renormalization. If needed, they would 
cause a finite shift in the Higgs-boson self-energies and thus would be relevant 
for precision calculations of M h . It is therefore important to clarify the necessity 
for such counterterms. 

In particular, if the symmetry is broken, the two-loop counterterms for the 
Higgs-potential have to be modified according to 

5 (2) V = 5^V sym + 5^V Qon _ susy , (1) 

which contains the usual multiplicative symmetric counterterms S^V sym and ad- 
ditional supersymmetry- restoring counterterms S^V non - susy . The renormalizabil- 
ity of the MSSM [30-33] guarantees that supersymmetry can be restored by a 
suitable choice of S^V non - susy . But since the structure of these extra counter- 
terms does not correspond to multiplicative renormalization, they have to be 
derived from the supersymmetric Slavnov- Taylor identities explicitly evaluated 
at the two-loop level. 

In the present paper we derive the relevant supersymmetric two-loop Slavnov- 
Taylor identities that are potentially affected by 5^V Don _ susy and determine the 
values of the supersymmetry-restoring counterterms by a direct calculation. Since 
we are aiming at contributions from the large Yukawa couplings, we restrict 
ourselves to the gauge-less limit with respect to the electroweak interaction, i.e. 
we take into account electroweak Yukawa interactions and the strong interaction 
in terms of supersymmetric QCD. 

Our explicit determination makes use of recent developments concerning di- 
mensional reduction [29, 34] (see also Refs. [35, 36] for complementary recent 
work). In Ref. [29], it has been shown that a slight reformulation of dimensional 
reduction, which does not affect its application to the calculation of Mh, avoids 
the mathematical inconsistency of Ref. [37]. Based on this reformulation, a dra- 
matic simplification for evaluating supersymmetric Slavnov- Taylor identities has 
been obtained. The evaluation of such identities constitutes the central part of 
our calculation of 5^V nan ^, KV . 
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The general strategy of our calculation follows the approach described in 
Refs. [25-27]. First we discuss general properties of the potential supersymmetry- 
restoring counterterms that are relevant for the calculation of the Higgs-boson 
mass. Then we derive suitable Slavnov- Taylor identities that are sensitive to 
these counterterms. Finally, the Slavnov- Taylor identities are evaluated using 
the method of Ref. [29]. 



2 Higgs potential and counterterms 

The two-loop evaluation of the Higgs-boson mass spectrum requires three types 
of counterterms: one-loop counterterms that renormalize divergent one-loop sub- 
diagrams, products of one-loop counterterms that partially renormalize two-loop 
diagrams, and genuine two-loop counterterms. As mentioned in the introduction, 
multiplicative renormalization is sufficient for the one-loop counterterms. Here 
we examine whether this applies also to the genuine two-loop counterterms. 

The genuine two-loop counterterms relevant for the evaluation of the Higgs- 
boson self-energies correspond to the MSSM Higgs potential V. At tree level, V 
is given by 

V = Quadratic "I" ^quartic; (^ a ) 
Quadratic = ™?l^l| 2 + ^|ft 2 | 2 + m\ (H\Hl - U\U\ + h.C.) , (2b) 

V = (l^il 2 - \n 2 \ 2 ) 2 + 1 \n\n 2 \ 2 (2c) 

in terms of the Higgs doublets Hi — H 1 + ("J), H2 = H 2 + ( ), their respective 
vacuum expectation values t>i )2 , and the SU(2) and U(l) gauge couplings g 2 ,i- 
The mass parameters m\ 2 are combinations of soft supersymmetry-breaking pa- 
rameters and the /z-parameter, and m\ is a soft supersymmetry-breaking param- 
eter. 

The symmetric counterterms used in the literature are generated by the mul- 
tiplicative renormalization transformation 

m^ — > m 2 + ^m? + ^m?, (3a) 

9i - 9i + 5^ 9l + 5^ 9i , (3b) 

Vi -> Vi + 5^v t + 5^v h (3c) 

H . _ (1 + S W ZBi + 5^Z Hi f/ 2 H h (3d) 

where the notation 8® refers to a renormalization constant of /-loop order. Thus, 
the genuine two-loop symmetric counterterms have the structure 



8 {2) V = 

u v sym 



1 z 



V, (4) 
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which is sufficient only if DRED preserves supersymmetry at the two-loop level 
in the gauge- less limit. 

If this is not the case, S^V has to be augmented by supersymmetry-restoring 
counterterms, as indicated in (JTJ). A priori, 5^ 2 'V non ^ BVSy contains quadratic and 
quartic terms, 

X(2)y _ r(2) ^quadratic , r(2) T^quartic /r\ 

u •'non-susy u * non— susy ^ u v non— susy \ u J 

Two simple properties of S^V noa - SUBy can immediately be derived. First, 
<5 (2) K on _ susy respects global SU(2) x U(l) gauge invariance. The reason is that 
DRED 1 preserves global gauge invariance as can easily be obtained e.g. from the 
gauge invariance of the regularized Lagrangian using the quantum action princi- 
ple [29]. As a consequence, 8^V non - susy depends on the Higgs doublets Hi 2 and 
vi t 2 only via the combinations Hip, similar to V itself. Second, the quadratic 
terms in 5^V non - susy are not necessary since the quadratic terms in the Higgs 
potential (|2bJ) are not restricted by supersymmetry in the first place, and no 
supersymmetry-breaking could modify their structure. Therefore, we have 

r(2)T / quadratic _ p / fi x 
u v non— susy w- V u / 



3 Relevant Slavnov- Taylor identities 

Any non-zero counterterm quartic in TCi, H.2 modifies both the Higgs-boson four- 
point and two-point vertex functions and can thus contribute to the Higgs- 
boson mass prediction. For example, the term |7ii| 4 = \H\\ A + 2v\\H\\ 2 + 
Av\{ReH\) 2 + . . . would contribute to the four-point function ^ HlH t HiH t and 
to the two-point functions involving H\. On the other hand, owing to van- 
ishing non-supersymmetric quartic counterterms would imply also the absence 
of non-supersymmetric counterterms for the self-energies. It is therefore suffi- 
cient to evaluate the quartic non-supersymmetric counterterms 5^ V^n-susy by 
considering their contribution to four-point functions. 

The starting point for the evaluation is the basic requirement that after 
renormalization, i.e. in particular after adding 5^V^^ sy , the supersymmet- 
ric Slavnov- Taylor identity of the MSSM [33] holds, 

S(T) = 0. (7) 

Here T denotes the renormalized effective action, the generating functional of the 
renormalized one-particle irreducible vertex functions. 

1 The version of DRED used in the literature is the one denned in Ref. [29] with anticom- 
muting 75. 
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The following derivative of the Slavnov- Taylor identity directly constrains the 
desired Higgs-boson four-point functions: 



5 5 S(T) 



SHLSeL 



(8) 



where fa denote any components of the MSSM Higgs bosons Hj, Hp, and H l k is 
the Higgsino partner of H k . e denotes the supersymmetry ghost. The index L 
denotes left-handed spinors, 



_5_ 

5e L 
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6 - 6 P 
5H[ L 5Hl L 



(9) 



with Pl = |(1 — 75). In the gauge-less limit, identity (JBJ) is equivalent to the 
following relation between renormalized vertex functions, 



+|r 

-E 

+ (!' 



a H l kL Y lj> .e I F<t>b<i>c<t>i 

r 



v , r 



3 

kh ±A i 



+ perm 



r ^ + r, jC . r 



H kL H i 



c + perm 



+ r , , jc- r _, ^rjc + perm 



ip .in 1 t~i th 

(paMcVitL Hi W 4>a<t>b<t>cVi C f-h Hi W° 

kL 1 kL i 

+V2P L f r (f>a(f>b(f>cHlkL -c. 



(10) 



Although this identity is quite involved, it is easy to see that it relates the de- 
sired quantity, the ^-interaction, which appears as the first term in the first 
line, to many other vertex functions. In the first sum, fa runs over all Higgs 
boson components H?, H(\ and in the second sum i,j take the values 1,2. The 
abbreviation "perm" denotes terms corresponding to all possible permutations 
of 4> a ,b,c- The conventions for the Slavnov- Taylor identity of the MSSM have 
been adapted to four-spinors as e.g. in Refs. [26,27]. The symbols Y^. and yj 
denote the sources of the BRS variations of fa and H?, respectively; x an d fo 
are the spinorial and constant components of the spurion superfield introduced 
to describe soft supersymmetry breaking [32]. 
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At the two- loop level, each term of the form Ta^b can be decomposed as 
follows, 



(0) ( -p(2),DR,ED 
B 



+ r 



(2),ct 
B 



+ r^rg* + ( r 



,(2),DRED 



+ r 



(2),ct^ r (0) 
B ' 



where r"' denotes the Z-loop order contribution. Thereby, the two-loop contribu- 
tion has been split into the regularized vertex functions including one-loop but 
excluding two- loop counterterms (denoted by the upper index "DRED"), and the 
corresponding genuine two-loop counterterms (index "ct"). 

Likewise, identity (JHJ) or (fTU|). evaluated at 2-loop order, can be rewritten in 
the form 



_ A (2),ct + A (2),DRED 



(12) 



where A( 2 )' ct is the contribution from genuine two-loop counterterms and 

A (2),DRED ig the 

remaining contribution, involving the regularized vertex func- 
tions and one-loop counterterms. In the gauge-less limit, the first contribution 
in (|T2*j) is given by 



A (2) 



ct 



£ 

<t>i 

£ 

hi 



r (o) r v 

H^Y^.tL <Pa<PbVc(Pi 



,(2),ct 



p(2),ct 



(0) 



i(2),ct 



r 



(0) 



_|_ p(2),ct pv 



,(0) 



+r 



(0) 



r 



(2),ct 



+ r 



(0) 



i(2),ct 



iC _ J. _ ?c + perm 



(13) 



All terms in (|1(J|) that do not re-appear in (JTHj) vanish since they involve Green 
functions of power-counting dimension > 5, for which there are neither classical 
nor counterterm contributions. Moreover, most of the terms in f!13|) vanish due 
to global gauge invariance or as a result of the gauge- less limit. The only term 
that can be non-zero is the first term for 6^ = 6 d = Hi. The result for T^) 

H kL Y 4>d tL 



is well known (see e.g. [26,27]), and we obtain 

A (2),ct = ^p^t 



(14) 



The second contribution in (fT2"|) has the same form as the right-hand side 
of (fTUJ) . but the products r^r^ have to be evaluated at the two-loop level without 
the genuine two-loop counterterms. A^ 2 ^ DRED can be written in the following 
simple form, 



A (2). 



DRED 



(5 5 5 -^ r (2),DRED^ 

*i(j) a 5(j) b 5(f) c 5(f) d 5e L 



(15) 



2— loop 
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where r( 2 - ) ' DRED denotes the regularized vertex functional including one-loop but 
excluding 2-loop counterterms. 

Hence, identity (JBJ) leads to 



Pt.T 



(2),ct 



1 

71 



DREDn 



■>cf><PdfcL J 2-loop 



(16) 



at the two-loop level in the gauge- less limit. This is an equation that directly 
determines the counterterm for any 4 -interaction. Since the multiplicative coun- 



terterm to rj? 
in the form 

5 4 



SfaSfaSfadtpa 



vanishes in the gauge-less limit, we can rewrite this equation 

6 5 S(T^ DRED ) 



d D xP L 8^V^Z\ 



non— susy 



1 

7! 



5<f) a 5<j) b 5<f) c 5<f) d 5e 1 



L J 2-loop 



(17) 



Obviously, by choosing all possible combinations of fa, this equation, together 
with its complex conjugate, completely determines the non-supersymmetric coun- 
terterms at the desired order. 



4 Evaluation of the Slavnov- Taylor identities 

The remaining task is to evaluate the right-hand side of (|17|) . Traditionally, we 
would have to evaluate all Green functions appearing in (|10|). in particular, also 
5-point functions at the two-loop level. We can simplify the evaluation by using 
the results of Ref. [29], according to which the right-hand side of (fTTj) is equal to 



A (2),DRED = _ %p 



([5(r(°))]-r DRED ) ( 



Pl, 

2— loop 



involving an insertion of the operator [S^T^)] at the two- loop level in the gauge- 
less limit. 

In Ref. [29], the operator [S(T^)] has been given for a general supersymmetric 
gauge theory without spontaneous breaking of gauge invariance and soft super- 
symmetry breaking. For the purpose of the present paper we have rederived this 
operator for the case of the MSSM, taking into account both symmetry breakings 
in the way defined in Ref. [33] . It turns out that those parts that can play a role 
in (J18|) are not modified compared to the case without symmetry breaking. All 
relevant terms in [S'(r^)] are four-fermion operators involving e and three other 
fermions (here either gluinos, Higgsinos, or quarks). 

Fig. Q shows the Feynman diagrams contributing in the gauge-less limit to 
(|18|) . i.e. to the vertex function with external <p a ,b,ci H l k , e and the insertion of the 
composite operator [S^I^ ))]. The insertion of [S^r^ ))] is marked by a cross, and 



S 



(a) (b) (c) 

Figure 1: Diagrams contributing to (JTSJl at the two-loop level in the gauge-less 
limit. The insertion of the operator [S'^ ))] is marked by a cross. Quarks, 
gluons and Higgs bosons are denoted by q, g, (ft; squarks, gluinos and Higgsinos 
are denoted by q, g, and H. The lines corresponding to the external Higgs bosons 
(fi a ,b,c have to be attached in all possible ways. 

the three Higgs fields have to be attached in all possible ways. In all diagrams, 
the basic fermion-loop topology is "Topology (c)" from Ref. [29]; i.e. the open 
fermion line involves e and quarks or Higgsinos, but no gauginos. 

It is useful to denote the chain of 7-matrices corresponding to the open fermion 
line as A, and the 7-chain corresponding to the closed fermion loop as B. Then 
the Feynman rules for the diagrams involve terms like the following ones, 

P Ll »B liX P L AP L , P L AP L Tr(P R B), ... (19) 

We do not need the detailed Feynman rules but only two properties that can be 
easily read off from the rules given in Ref. [29] combined with the explicit factors 
Pl appearing in (fTKj) : 

1. a diagram vanishes if the number of 7^-matrices in AB is odd. 

2. a diagram vanishes if the number of 7 M -matrices in B is smaller than four. 

After carrying out the two loop integrals in diagrams of the form Fig. [TJi,b, the 
only covariants that can appear in the 7-chain B are tf abc , where p a ,b,c are the 
incoming momenta of the Higgs fields <p a ,b,c- Hence the 7-chain B can be simplified 
to terms involving at most three factors of and no other 7^-matrices. As a 
result of property 2 above, these diagrams all vanish. 

Similarly, it can be easily seen that diagrams of the form Fig. vanish if 
at most two of the three Higgs fields are attached to the fermion loop. The 
only remaining case is the one shown in Fig. El where all three Higgs fields are 
attached to the fermion loop. This case is now discussed in more detail. Due to 
the structure ^ + rrii of the numerators of the fermion propagators, the 7-chain B 
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Figure 2: Diagram of the form shown in Fig. where all Higgs-boson lines are 
attached to the fermion loop. In this instance the 7-chains A and B are products 
of one and five fermion propagators, respectively. 

in the diagram of Fig- EI is a sum of terms involving between zero and five factors 
of d v while the terms of the 7-chain A involve either zero or one such factor. 

After carrying out the fermion-loop integral, B can contain up to four covari- 
ants </j G {]f> abc1 ft}, where k is the second loop momentum. Hence in those terms 
of B that contain five factors of (j v actually two of the ^ must be equal, and 
the product of these five ^ can be reduced to products of only three different tf v 
Hence such terms of B contribute zero to the diagrams. 

Those terms of B that contain four factors of ^ can contribute only together 
with terms of A containing no such factor, due to property 1 above. In such a 
situation, the fc-dependence of the second loop integral has the following form, 

d D k ft f(k 2 , kp a , kp b , kp c ) 



21 ' 



(2tt) d \(k + p a +Pb + Pc) 2 - m 2 ] [k 2 - m 

where the scalar function / results from the fermion-loop integral. As a result of 
the loop integral (J2*U|) . the factor ft within B is replaced by a linear combination 
of f' a b . In this way, each term within B can be simplified to a term involving at 
most two ^-factors. 

Taken together, there are no terms in the 7-chain B that can possibly give a 
contribution of diagrams of the form Fig. Therefore, all diagrams in Fig. Q 
vanish, and thus we obtain 

A (2),DRED = 0; (21) 

and, with the help of (|15|) and (J17)) . we find the relation 

5 {2) VZ:^ sy = , (22) 
which excludes the presence of non-symmetric counterterms. 
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5 Conclusions 



We have considered the Yukawa-enhanced two-loop contributions to Higgs-boson 
masses in the MSSM, i.e. contributions of 0{a t ^a s ), 0(af b ), 0(a t ab). In the 
literature these contributions have been evaluated under the assumption that 
DRED is a supersymmetric regularization up to the required order. In the present 
article we have verified this assumption. 

The question of necessity for supersymmetry-restoring counterterms was stud- 
ied utilizing the algebraic method of supersymmetric Slavnov- Taylor identities. 
A Slavnov- Taylor identity has been found that unambiguously determines the 
potentially necessary non-symmetric counterterms. It involves up to two-loop 5- 
point functions, but it could be explicitly evaluated using the method of Ref. [29]. 
The identity turned out to be valid at the regularized level, which means that su- 
persymmetry is preserved by DRED, multiplicative renormalization is sufficient, 
and therefore no extra supersymmetry-restoring counterterms are necessary. 

The restriction to the Yukawa-enhanced contributions and the gauge-less limit 
is of crucial importance. If the gauge-less limit is relaxed, more terms of the 
operator [S(T^)) can contribute and can lead to more complicated Feynman 
diagrams. For example, there are diagrams with a topology like in Fig. \I]p, but 
with a vector boson instead of the virtual squark (and e.g. the if-g-g-vertex 
replaced by a H-H-V- vertex) . The 7-strings A and B in such diagrams have 
one 7-matrix more, and the arguments of section 0] cannot be applied. 

It is, therefore, not straightforward to extrapolate the results of the present 
paper to the full electroweak two- loop contributions and/or to the three- loop 
level. The future experimental accuracy of Higgs-boson mass measurements, 
however, requires to bring the full two-loop and leading three-loop contributions 
under control. A dedicated study of the supersymmetry-properties of DRED for 
these more complicated cases will be an indispensable step in this program. 
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